This project seeks to estimate sport fish harvest and releases of rockfish in Alaska waters by improving on the Howard et al. (2020) methods and expand the time series back to 1977 when the statewide harvest survey (SWHS) was first implemented. This is essentially a Bayesian version of the Howard methods that allows for more appropriate and defensible sharing of information between areas, handles missing data in a more appropriate manor, accurately propagates uncertainty throughout the estimation procedure and replaces the Howard decision tree approach to low sample sizes with a hierarchical model. The methods and results for generating harvest estimates are generally consistent between the Bayesian model and the Howard methods. Harvest estimates are consistent with Howard estimates during contemporary times, but may differ based on more appropriate weighting of SWHS and logbook data, including estimating and correcting bias in the SWHS data.

Table 1. Summary of key improvements in reconstructiing sport fish removals of rockfish using the Bayesian model as compared to the Howard et al. (2020) methods.
Issue Howard Bayes
Time series 1999 - present 1977 - present
Bias in SWHS Not explicitly dealt with. Relies on logbook data and ratios of guided/unguided from SWHS data to estimate unguided releases and harvests. Explicitly estimates bias in SWHS harvest and release estimates based on logbook data.
Species composition of releases Assumes that species composition of releases is equal to that of the harvest, which is not evident in the logbook data. Recognizes different release probabilities by species / species assemblage and estimates it from logbook data and bias corrected SWHS data
Sample size limitations Uses sample size threshholds such that when areas fall below those threshholds values are borrowed from nearby areas. Uses a hierarchichacal modelling approach that shares information between areas in the same region. Thus all data is used, even with small sample sizes. This is a more sound method that avoids assumptions and uses all of the data.
Error propogation Error is propogated when variance estimates are available, but there is uncertainty associated with borrowing values from nearby areas, or the assumption of species compositions being identical in harvest and releases, are not dealt with. By breaking the assumption that species composition is equal between harvests and releases, uncertainty in the release estimates is more reflective of the fishery. Furthermore, the hyerarchichal approach more accurately captures uncertainy within and between areas within a region.
**Figure 1.**- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.

Figure 1.- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.


**Figure 2.**- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).

Figure 2.- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).


Results

**Figure X.**- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Figure X.- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Estimate comparison

Since previous estimates of rockfish harvest have been produced these first 3 graphs will be used to show how the modeled estimates compare to the estimates produced earlier. For total rockfish the estimates are in general agreement although differences are noted. These estimates should be more reliable because they include both SWHS and guide logbook data, handle variance more appropriately, use hierarchical distributions when data is missing, directly consider observation error and are produced using reproducible research.

**Figure 3.**- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 3.**- Total rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.


Notes from Adam: When looking at only black rockfish the most significant differences are for the Prince William Sound Inside area. I did not spend a great deal of time tracking this down although it looks like the previous version used bad values for \(P_{(black)ayu}\) for at least unguided anglers. For the moment I would ignore the results for BSIA and SOKO2SAP. I think it is possible to give approximate values for these areas but it will require a little more coding which I have yet to do.

**Figure 4.**- Black rockfish harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 4.- Black rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.


And black rockfish releases…

**Figure 5.**- Black rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 5.- Black rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.





**Figure 6.**- Yellow rockfish harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 6.- Yellow rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 7.**- Yellow rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 7.- Yellow rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.





**Figure 8.**- DSR rockfish (excluding yelloweye) harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 8.- DSR rockfish (excluding yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 9.**- DSR rockfish releases (including yelloweye) 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 9.- DSR rockfish releases (including yelloweye) 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 11.**- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 11.- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 12.**- Slope rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 12.- Slope rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

**Figure 13.**- *Other* rockfish (non-pelagic and non-yelloweye) harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 13.- Other rockfish (non-pelagic and non-yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 14.**- *Other* rockfish (non-pelagic and non-yelloweye) releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 14.- Other rockfish (non-pelagic and non-yelloweye) releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



Total Biomass Removal Estimates

**Figure 15.**- Black rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 15.- Black rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.



**Figure 16.**- Yellow rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 16.- Yellow rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

**Figure 17.**- Pelagic rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 17.- Pelagic rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


**Figure 18.**- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 18.- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


**Figure 19.**- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 19.- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


Model fit

Logbook residuals

**Figure 20.**- Residuals from logbook harvests.

Figure 20.- Residuals from logbook harvests.


SWHS residuals

**Figure 21.**- Residuals from SWHS harvests.

Figure 21.- Residuals from SWHS harvests.



**Figure 22.**- Residual of SWHS releases.

Figure 22.- Residual of SWHS releases.

Parameter estimates

P(Charter)

These histograms show the posterior distribution of the mean percent of rockfish harvested by the charter fleet.

**Figure 23.**- Mean percent of harvest by charter anglers.

Figure 23.- Mean percent of harvest by charter anglers.


When considered annually we see the percent of rockfish harvested by the charter fleet follows our data fairly well although the model smooths out the changes and we just do not have much information about this ratio. Prior to 2011 the percent charter is confounded with SWHS bias and should be mostly discounted.

**Figure 24.**- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

Figure 24.- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

P(Harvest)

These plots show the fitted logistic line to the proportion of caught rockfish that are harvested. These estimates are used for hindcasting catch estimates based on the harvest data in early years when catch estimates are unavailable.


**Figure 25.**- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 25.- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.


**Figure 26.**- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 26.- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.


**Figure 27.**- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.

Figure 27.- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.


**Figure 28.**- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 28.- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.


**Figure 29.**- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 29.- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

SWHS bias

Figure 23 shows the mean estimate for SWHS bias in harvests and releases. Cook Inlet, North Gulf Coast and North Southeast Inside all look pretty good while most other areas have substantial bias. Prince William Sound Inside has the largest bias. Bias in release estimates is substantial and whereas the SWHS appears to underestimate harvests, it appears to greatly overestimates releases by a factor of 2 or more in most areas as derived from logbook reported releases.

**Figure 30.**- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS *underestimates* the true value and bias > 1 indicates the survey *overestimates* the true value.

Figure 30.- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.


Our estimates of SWHS harvest bias track observations fairly well when he have guided harvest estimates. The estimates of release bias in the SWHS data track observed patterns to an extent, but appear to smooth these more volatile disagreements with the logbook data. Adam postulated in his initial start on this that some of this could be the result of the estimates of the proportion guided. This value was not modelled with a trend and thus applies a constant estimate when hindcasting. Data on these relationships could greatly improve this model.

**Figure 31.**- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS *underestimates* the true value and bias > 1 indicates the survey *overestimates* the true value.

Figure 31.- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.

P(pelagic)

We model the percentage of pelagic rockfish in the harvest because we have the information for charter anglers (via logbooks) starting in 1998. Other than looking at the model estimates you can use Figure 25 to compare the two data streams for pelagic rockfish harvest. In general they are in agreement with major exceptions in Price William Sound inside, Prince William Sound outside (early in the time series) and South Southeast inside.

**Figure 32.**- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

Figure 32.- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

P(black|pelagic)

Note that in Southeast Alaska we only have composition data starting in 2006. Tania dug up old SE data, but it did not provide any useful data for species apportionment. For the most part, P(black|pelagic) is relatively constant across areas, with the exception of Cook Inlet and NSEI in Southeast AK. It may be worth discussing whether the shifts in those areas is a result of improved or changing species identification rather than actual shift in the species composition of the catch.

**Figure 33.**- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

Figure 33.- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

P(yelloweye|non-pelagic / yelloweye|DSR)

**Figure 34.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

Figure 34.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

P(DSR|non-pelagic)

**Figure 35.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

Figure 35.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

P(slope|non-pelagic)

**Figure 36.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

Figure 36.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.



P(slope|non-pelagic & non-yellowye) For release estimates

## NULL



Weight Fits

**Figure 37.**- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 37.- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 38.**- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 38.- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 39.**- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 39.- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 40.**- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 40.- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 41.**- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 41.- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


Summary of unconverged parameters:

The proportion of unconverged species composition proportional estimates that are near the bounds of 0 or 1 is 0.3515152

The proportion of unconverged true parameters (defined as Rhat > 1.01) that are also in excess of an Rhat of 1.1 is 0.1 while 0.9 are less than 1.1.

Table Xa. Species proportion derived quantities that are not with 0.05 of the 0 and 1 bounds and are not converged relative to an Rhat of 1.01.

variable

area_n

year

user

mean

median

lower

upper

Rhat

sd

skew

area

p_pelagic

11

12

2

0.395

0.388

0.182

0.651

1.010

0.122

0.057

CSEO

p_pelagic

11

15

1

0.434

0.429

0.211

0.679

1.020

0.124

0.040

CSEO

p_pelagic

11

18

1

0.477

0.474

0.247

0.716

1.011

0.122

0.025

CSEO

p_pelagic

13

1

1

0.310

0.296

0.104

0.600

1.086

0.131

0.107

NSEI

p_pelagic

13

1

2

0.122

0.107

0.031

0.310

1.075

0.071

0.210

NSEI

p_yellow

13

1

1

0.421

0.384

0.135

0.921

1.012

0.197

0.188

NSEI

p_yellow

13

1

2

0.385

0.338

0.119

0.907

1.013

0.197

0.239

NSEI

p_pelagic

13

10

1

0.330

0.316

0.120

0.604

1.050

0.126

0.111

NSEI

p_pelagic

13

10

2

0.135

0.119

0.041

0.312

1.060

0.071

0.224

NSEI

p_yellow

13

10

1

0.418

0.386

0.143

0.916

1.010

0.191

0.168

NSEI

p_yellow

13

10

2

0.382

0.342

0.118

0.904

1.016

0.194

0.206

NSEI

p_pelagic

13

11

1

0.332

0.318

0.128

0.602

1.052

0.125

0.112

NSEI

p_pelagic

13

11

2

0.135

0.120

0.041

0.315

1.054

0.072

0.209

NSEI

p_yellow

13

11

1

0.416

0.377

0.132

0.915

1.010

0.195

0.200

NSEI

p_yellow

13

11

2

0.389

0.347

0.121

0.897

1.012

0.194

0.216

NSEI

p_pelagic

13

12

1

0.341

0.330

0.133

0.609

1.035

0.126

0.087

NSEI

p_pelagic

13

12

2

0.139

0.125

0.042

0.321

1.041

0.072

0.193

NSEI

p_yellow

13

12

2

0.385

0.340

0.120

0.898

1.012

0.196

0.230

NSEI

p_pelagic

13

13

1

0.342

0.330

0.135

0.619

1.033

0.124

0.097

NSEI

p_pelagic

13

13

2

0.140

0.125

0.046

0.310

1.045

0.071

0.211

NSEI

p_yellow

13

13

1

0.423

0.388

0.134

0.916

1.012

0.196

0.179

NSEI

p_yellow

13

13

2

0.385

0.343

0.122

0.898

1.011

0.193

0.217

NSEI

p_pelagic

13

14

1

0.327

0.311

0.133

0.599

1.037

0.120

0.133

NSEI

p_pelagic

13

14

2

0.139

0.127

0.045

0.307

1.029

0.068

0.175

NSEI

p_pelagic

13

15

1

0.302

0.288

0.122

0.570

1.047

0.116

0.121

NSEI

p_pelagic

13

15

2

0.132

0.119

0.045

0.287

1.036

0.063

0.205

NSEI

p_yellow

13

15

1

0.419

0.383

0.138

0.918

1.011

0.195

0.185

NSEI

p_yellow

13

15

2

0.379

0.339

0.118

0.907

1.012

0.194

0.206

NSEI

p_pelagic

13

16

1

0.340

0.325

0.141

0.588

1.023

0.119

0.126

NSEI

p_pelagic

13

16

2

0.145

0.132

0.047

0.310

1.027

0.069

0.189

NSEI

p_yellow

13

16

2

0.387

0.347

0.121

0.896

1.010

0.195

0.205

NSEI

p_pelagic

13

17

1

0.352

0.340

0.161

0.603

1.018

0.116

0.103

NSEI

p_pelagic

13

17

2

0.150

0.138

0.052

0.317

1.030

0.070

0.172

NSEI

p_yellow

13

17

2

0.384

0.345

0.121

0.894

1.012

0.193

0.202

NSEI

p_pelagic

13

18

1

0.361

0.350

0.161

0.608

1.017

0.116

0.095

NSEI

p_pelagic

13

18

2

0.154

0.143

0.055

0.323

1.033

0.070

0.156

NSEI

p_yellow

13

18

2

0.381

0.341

0.119

0.899

1.013

0.194

0.207

NSEI

p_pelagic

13

19

2

0.159

0.148

0.055

0.334

1.012

0.070

0.158

NSEI

p_yellow

13

19

1

0.417

0.381

0.136

0.912

1.013

0.193

0.187

NSEI

p_yellow

13

19

2

0.380

0.338

0.119

0.895

1.011

0.192

0.218

NSEI

p_pelagic

13

2

1

0.311

0.291

0.108

0.606

1.079

0.134

0.149

NSEI

p_pelagic

13

2

2

0.122

0.106

0.032

0.295

1.069

0.069

0.232

NSEI

p_yellow

13

2

1

0.421

0.386

0.142

0.913

1.010

0.193

0.182

NSEI

p_yellow

13

2

2

0.386

0.343

0.120

0.914

1.012

0.195

0.220

NSEI

p_pelagic

13

20

1

0.373

0.366

0.177

0.618

1.018

0.116

0.060

NSEI

p_pelagic

13

20

2

0.167

0.155

0.058

0.339

1.010

0.073

0.163

NSEI

p_yellow

13

20

2

0.383

0.340

0.124

0.898

1.011

0.195

0.221

NSEI

p_yellow

13

21

2

0.383

0.341

0.119

0.889

1.012

0.195

0.216

NSEI

p_yellow

13

22

1

0.418

0.382

0.138

0.918

1.011

0.194

0.186

NSEI

p_yellow

13

22

2

0.384

0.345

0.114

0.891

1.012

0.192

0.203

NSEI

p_yellow

13

23

1

0.431

0.401

0.145

0.908

1.011

0.192

0.156

NSEI

p_yellow

13

23

2

0.382

0.341

0.117

0.886

1.010

0.191

0.215

NSEI

p_yellow

13

24

1

0.434

0.401

0.151

0.906

1.014

0.191

0.172

NSEI

p_yellow

13

25

1

0.437

0.403

0.148

0.894

1.010

0.189

0.180

NSEI

p_yellow

13

27

2

0.368

0.335

0.119

0.869

1.012

0.183

0.181

NSEI

p_yellow

13

28

1

0.425

0.406

0.144

0.847

1.012

0.180

0.106

NSEI

p_pelagic

13

3

1

0.314

0.300

0.111

0.597

1.066

0.129

0.109

NSEI

p_pelagic

13

3

2

0.124

0.107

0.032

0.312

1.066

0.073

0.232

NSEI

p_yellow

13

3

1

0.419

0.381

0.141

0.915

1.011

0.195

0.195

NSEI

p_yellow

13

3

2

0.383

0.339

0.114

0.900

1.012

0.197

0.223

NSEI

p_pelagic

13

4

1

0.315

0.298

0.111

0.598

1.080

0.130

0.131

NSEI

p_pelagic

13

4

2

0.126

0.109

0.035

0.301

1.071

0.070

0.243

NSEI

p_yellow

13

4

2

0.390

0.350

0.117

0.904

1.013

0.196

0.204

NSEI

p_pelagic

13

5

1

0.318

0.302

0.112

0.605

1.072

0.128

0.125

NSEI

p_pelagic

13

5

2

0.126

0.111

0.034

0.304

1.076

0.072

0.208

NSEI

p_yellow

13

5

1

0.420

0.384

0.140

0.914

1.011

0.194

0.186

NSEI

p_yellow

13

5

2

0.382

0.341

0.121

0.897

1.012

0.195

0.210

NSEI

p_pelagic

13

6

1

0.319

0.305

0.112

0.606

1.062

0.127

0.110

NSEI

p_pelagic

13

6

2

0.128

0.113

0.037

0.306

1.059

0.070

0.215

NSEI

p_yellow

13

6

1

0.420

0.384

0.135

0.918

1.017

0.195

0.184

NSEI

p_yellow

13

6

2

0.387

0.349

0.121

0.915

1.015

0.195

0.195

NSEI

p_pelagic

13

7

1

0.321

0.307

0.118

0.609

1.071

0.130

0.107

NSEI

p_pelagic

13

7

2

0.129

0.113

0.036

0.316

1.056

0.072

0.223

NSEI

p_yellow

13

7

2

0.381

0.341

0.110

0.894

1.015

0.194

0.206

NSEI

p_pelagic

13

8

1

0.324

0.311

0.123

0.599

1.062

0.127

0.103

NSEI

p_pelagic

13

8

2

0.129

0.114

0.037

0.296

1.056

0.069

0.217

NSEI

p_yellow

13

8

1

0.421

0.387

0.144

0.907

1.012

0.193

0.176

NSEI

p_yellow

13

8

2

0.381

0.338

0.121

0.898

1.013

0.194

0.222

NSEI

p_pelagic

13

9

1

0.324

0.307

0.120

0.601

1.053

0.128

0.133

NSEI

p_pelagic

13

9

2

0.132

0.116

0.038

0.307

1.050

0.071

0.226

NSEI

p_yellow

13

9

1

0.420

0.385

0.135

0.924

1.015

0.195

0.180

NSEI

p_yellow

13

9

2

0.387

0.346

0.123

0.898

1.016

0.193

0.212

NSEI

p_pelagic

16

11

1

0.316

0.305

0.128

0.559

1.014

0.112

0.098

SSEO

p_pelagic

16

12

1

0.310

0.300

0.123

0.562

1.014

0.113

0.088

SSEO

p_pelagic

16

13

2

0.216

0.201

0.077

0.424

1.012

0.092

0.164

SSEO

p_pelagic

16

14

1

0.295

0.284

0.116

0.540

1.020

0.107

0.103

SSEO

p_pelagic

16

15

1

0.299

0.290

0.122

0.529

1.013

0.107

0.084

SSEO

p_pelagic

16

16

1

0.295

0.284

0.121

0.524

1.012

0.106

0.104

SSEO

p_pelagic

16

18

2

0.221

0.206

0.076

0.446

1.011

0.095

0.157

SSEO

p_pelagic

16

20

1

0.315

0.304

0.141

0.545

1.013

0.105

0.105

SSEO

p_pelagic

16

4

1

0.313

0.305

0.116

0.559

1.013

0.114

0.070

SSEO

p_pelagic

16

4

2

0.218

0.206

0.076

0.440

1.010

0.092

0.130

SSEO

p_pelagic

16

5

1

0.313

0.302

0.123

0.550

1.011

0.112

0.098

SSEO

p_pelagic

16

6

1

0.310

0.300

0.123

0.560

1.011

0.113

0.088

SSEO

p_pelagic

16

8

1

0.311

0.302

0.122

0.559

1.010

0.113

0.079

SSEO

p_pelagic

16

8

2

0.216

0.203

0.073

0.430

1.012

0.092

0.141

SSEO

p_pelagic

3

1

1

0.542

0.547

0.272

0.794

1.012

0.136

-0.037

PWSI

p_pelagic

3

11

1

0.546

0.548

0.283

0.782

1.010

0.132

-0.015

PWSI

p_pelagic

3

12

1

0.542

0.550

0.276

0.780

1.022

0.133

-0.060

PWSI

p_pelagic

3

14

1

0.548

0.557

0.264

0.784

1.011

0.133

-0.067

PWSI

p_pelagic

3

14

2

0.185

0.171

0.061

0.388

1.013

0.084

0.167

PWSI

p_pelagic

3

15

1

0.542

0.549

0.276

0.780

1.011

0.131

-0.054

PWSI

p_pelagic

3

2

2

0.182

0.168

0.060

0.377

1.012

0.082

0.170

PWSI

p_pelagic

3

3

1

0.542

0.549

0.274

0.784

1.010

0.133

-0.053

PWSI

p_pelagic

3

3

2

0.180

0.166

0.055

0.380

1.016

0.084

0.167

PWSI

p_pelagic

3

8

1

0.543

0.550

0.273

0.780

1.017

0.134

-0.052

PWSI

p_black

8

40

1

0.432

0.409

0.339

0.930

1.015

0.115

0.200

afognak

Table Xb. True parameters with an Rhat in excess of 1.1

variable

area_n

year/species

mean

median

lower

upper

Rhat

sd

skew

year

area

beta2_pH

10

1

1.324

0.753

0.274

6.786

1.205

1.592

0.359

northeast

beta0_pH

11

1

0.090

0.126

-0.638

0.700

1.247

0.347

-0.104

CSEO

beta1_pH

11

1

3.303

3.218

2.367

4.574

1.271

0.556

0.153

CSEO

beta2_pH

11

1

0.429

0.370

0.178

1.077

1.107

0.257

0.230

CSEO

beta2_yellow

11

-6.489

-4.967

-20.216

-0.260

1.151

5.314

-0.286

CSEO

beta0_pelagic

13

-1.002

-0.973

-2.115

-0.127

1.267

0.529

-0.055

NSEI

beta1_pelagic

13

3.409

3.137

1.690

6.379

1.167

1.195

0.228

NSEI

beta2_pelagic

13

0.154

0.140

0.064

0.343

1.199

0.077

0.182

NSEI

beta0_pelagic

16

-0.821

-0.789

-1.460

-0.442

1.137

0.238

-0.135

SSEO

beta1_pelagic

16

1.954

1.907

1.412

2.929

1.122

0.353

0.133

SSEO

mu_beta2_pH

2

1

1.043

0.501

0.034

5.877

1.266

1.609

0.337

NG

tau_beta2_pH

2

3

0.256

0.109

0.042

1.231

1.123

0.773

0.190

NG

beta2_pH

6

1

0.965

0.352

0.103

6.458

1.149

1.627

0.377

SOKO2SAP

beta2_pH

8

1

0.827

0.288

0.111

6.818

1.350

1.645

0.328

afognak

beta2_pH

9

1

1.234

0.494

0.152

7.310

1.333

1.849

0.400

eastside

beta2_pH

9

3

0.188

0.159

0.081

0.429

1.135

0.165

0.175

eastside

Table Xc. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.01.

Statistic

non.BSAI.SOKO.Harv

BSAI.SOKO.Harv

non.BSAI.SOKO.Rel

BSAI.SOKO.Rel

n unconverged at Rhat = 1.01

399.00

246.00

947.00

335.00

Prop. of total

0.62

0.38

0.74

0.26

Prop. median less than 1

0.28

0.20

0.33

0.43

Prop. median less than 10

0.41

0.57

0.52

0.77

Prop. median less than 50

0.50

0.79

0.65

0.83

Median of the median values

52.03

7.55

6.93

1.71

Mean of the median values

492.27

352.15

314.41

201.04

Mean skew of posteriors

0.24

0.29

0.27

0.25

Table Xd. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.1.

Statistic

non.BSAI.SOKO.Harv

BSAI.SOKO.Harv

non.BSAI.SOKO.Rel

BSAI.SOKO.Rel

n unconverged at Rhat = 1.1

7.00

8.00

42.00

34.00

Prop. of total

0.47

0.53

0.55

0.45

Prop. median less than 1

1.00

0.62

0.88

0.65

Prop. median less than 10

1.00

1.00

0.95

1.00

Prop. median less than 50

1.00

1.00

0.98

1.00

Median of the median values

0.00

0.66

0.00

0.21

Mean of the median values

0.06

1.09

4.31

0.84

Mean skew of posteriors

0.18

0.20

0.15

0.15

Table Xe. All parameters and derived quantities that are not converged relative to an Rhat of 1.1.

variable

Freq

beta0_pelagic

2

beta0_pH

1

beta1_pelagic

2

beta1_pH

1

beta2_pelagic

1

beta2_pH

5

beta2_yellow

1

eps_pH

3

H_ayg

1

Hb_ay

1

Hb_ayg

1

Ho_ay

3

Ho_ayu

8

Hp_ay

1

Hp_ayg

1

Hs_ayu

1

Hy_ay

1

Hy_ayu

2

mu_beta2_pH

1

p_slope

1

R_ayg

2

R_ayu

1

Rb_ayg

2

Rd_ay

1

Rd_ayg

1

Rd_ayu

2

re_pH

3

Ro_ay

8

Ro_ayg

3

Ro_ayu

27

Rp_ayg

2

Rs_ay

1

Rs_ayg

1

Rs_ayu

2

Ry_ay

11

Ry_ayg

9

Ry_ayu

22

tau_beta2_pH

1

Parameter estimates:

Summary Table of Parameter Estimates
Parameter mean sd Lower_CI Median Upper_CI
mu_bc_H[1] -0.021 0.069 -0.147 -0.024 0.122
mu_bc_H[2] 0.076 0.054 -0.022 0.073 0.190
mu_bc_H[3] -0.268 0.068 -0.396 -0.269 -0.129
mu_bc_H[4] -0.923 0.180 -1.280 -0.916 -0.569
mu_bc_H[5] 1.345 1.233 0.101 1.061 4.201
mu_bc_H[6] -2.367 0.320 -2.968 -2.366 -1.741
mu_bc_H[7] -0.456 0.103 -0.665 -0.453 -0.262
mu_bc_H[8] 0.269 0.337 -0.311 0.245 0.995
mu_bc_H[9] -0.248 0.126 -0.493 -0.249 0.009
mu_bc_H[10] 0.060 0.067 -0.066 0.060 0.195
mu_bc_H[11] -0.122 0.038 -0.197 -0.122 -0.049
mu_bc_H[12] -0.249 0.105 -0.473 -0.243 -0.050
mu_bc_H[13] -0.153 0.073 -0.287 -0.153 -0.010
mu_bc_H[14] -0.303 0.089 -0.477 -0.301 -0.135
mu_bc_H[15] -0.330 0.046 -0.419 -0.330 -0.239
mu_bc_H[16] -0.209 0.367 -0.836 -0.229 0.592
mu_bc_R[1] 1.032 0.121 0.793 1.032 1.276
mu_bc_R[2] 1.011 0.077 0.853 1.015 1.152
mu_bc_R[3] 0.229 0.091 0.044 0.231 0.405
mu_bc_R[4] 0.666 0.207 0.225 0.673 1.046
mu_bc_R[5] 1.674 0.365 0.951 1.671 2.409
mu_bc_R[6] -1.425 0.351 -2.107 -1.426 -0.753
mu_bc_R[7] 0.517 0.167 0.191 0.522 0.841
mu_bc_R[8] 0.638 0.150 0.344 0.634 0.941
mu_bc_R[9] 0.377 0.169 0.014 0.382 0.695
mu_bc_R[10] 0.900 0.106 0.690 0.901 1.105
mu_bc_R[11] 1.112 0.078 0.963 1.109 1.270
mu_bc_R[12] 0.894 0.168 0.563 0.898 1.223
mu_bc_R[13] 1.059 0.091 0.875 1.060 1.233
mu_bc_R[14] 0.842 0.126 0.598 0.840 1.090
mu_bc_R[15] 0.728 0.095 0.541 0.727 0.913
mu_bc_R[16] 1.302 0.106 1.088 1.304 1.506
tau_pH[1] 2.445 0.179 2.100 2.439 2.811
tau_pH[2] 1.568 0.146 1.297 1.563 1.876
tau_pH[3] 2.283 0.304 1.731 2.259 2.927
tau_pH[4] 1.050 0.143 0.797 1.048 1.353
tau_pH[5] 1.028 0.154 0.756 1.017 1.360
beta0_pH[1,1] 0.546 0.176 0.195 0.551 0.884
beta0_pH[2,1] 1.428 0.165 1.097 1.429 1.728
beta0_pH[3,1] 1.401 0.215 0.906 1.424 1.748
beta0_pH[4,1] 1.554 0.236 1.038 1.564 1.970
beta0_pH[5,1] -0.843 0.400 -1.814 -0.802 -0.186
beta0_pH[6,1] -0.635 0.523 -2.083 -0.508 0.029
beta0_pH[7,1] 0.752 0.011 0.731 0.752 0.773
beta0_pH[8,1] -0.773 0.387 -1.733 -0.720 -0.153
beta0_pH[9,1] -0.336 0.285 -0.990 -0.311 0.148
beta0_pH[10,1] 0.261 0.197 -0.148 0.265 0.616
beta0_pH[11,1] 0.090 0.347 -0.638 0.126 0.700
beta0_pH[12,1] 0.645 0.193 0.255 0.654 0.999
beta0_pH[13,1] -0.646 0.140 -0.920 -0.651 -0.372
beta0_pH[14,1] 0.067 0.180 -0.287 0.067 0.413
beta0_pH[15,1] 0.070 0.188 -0.349 0.083 0.406
beta0_pH[16,1] -0.963 0.222 -1.405 -0.954 -0.548
beta0_pH[1,2] 2.567 0.084 2.404 2.567 2.730
beta0_pH[2,2] 3.078 0.022 3.037 3.079 3.121
beta0_pH[3,2] 3.625 0.041 3.546 3.626 3.704
beta0_pH[4,2] 2.936 0.019 2.899 2.937 2.974
beta0_pH[5,2] 4.501 1.419 2.681 4.209 8.002
beta0_pH[6,2] 3.326 0.160 3.020 3.318 3.646
beta0_pH[7,2] 2.022 0.050 1.927 2.023 2.120
beta0_pH[8,2] 2.823 0.064 2.700 2.821 2.948
beta0_pH[9,2] 3.963 0.204 3.574 3.961 4.377
beta0_pH[10,2] 4.007 0.131 3.762 4.004 4.277
beta0_pH[11,2] 1.888 0.009 1.871 1.888 1.906
beta0_pH[12,2] 1.802 0.068 1.674 1.802 1.940
beta0_pH[13,2] 2.498 0.026 2.447 2.498 2.549
beta0_pH[14,2] 1.805 0.025 1.757 1.805 1.854
beta0_pH[15,2] 2.550 0.017 2.518 2.550 2.584
beta0_pH[16,2] 2.777 0.025 2.729 2.777 2.826
beta0_pH[1,3] 0.686 0.602 -1.013 0.849 1.365
beta0_pH[2,3] 2.126 0.018 2.091 2.126 2.161
beta0_pH[3,3] 2.690 0.023 2.644 2.690 2.736
beta0_pH[4,3] 2.980 0.033 2.916 2.980 3.046
beta0_pH[5,3] 3.826 2.680 1.028 3.134 11.439
beta0_pH[6,3] 1.256 0.322 0.691 1.232 1.954
beta0_pH[7,3] 0.856 0.049 0.760 0.857 0.952
beta0_pH[8,3] 0.283 0.123 0.011 0.292 0.502
beta0_pH[9,3] -2.859 1.091 -5.253 -2.747 -1.025
beta0_pH[10,3] -0.474 0.756 -2.252 -0.358 0.682
beta0_pH[11,3] -0.429 11.328 -21.744 -0.638 22.052
beta0_pH[12,3] -0.522 11.486 -23.988 -0.677 21.749
beta0_pH[13,3] -0.445 11.371 -22.672 -0.543 21.930
beta0_pH[14,3] -0.376 11.381 -22.930 -0.394 22.345
beta0_pH[15,3] -0.506 11.397 -23.836 -0.480 21.854
beta0_pH[16,3] -0.339 11.614 -22.672 -0.438 22.523
beta0_pH[11,4] 1.383 0.980 -1.600 1.551 2.821
beta0_pH[12,4] 1.366 1.074 -1.239 1.393 3.318
beta0_pH[13,4] 1.859 0.783 0.072 1.928 3.151
beta0_pH[14,4] 0.651 0.678 -1.011 0.733 1.748
beta0_pH[15,4] 0.580 0.964 -1.937 0.928 1.607
beta0_pH[16,4] 1.168 0.872 -0.944 1.284 2.613
beta0_pH[11,5] 1.289 0.113 1.096 1.279 1.535
beta0_pH[12,5] -2.311 0.430 -3.275 -2.275 -1.596
beta0_pH[13,5] 0.762 0.054 0.664 0.759 0.883
beta0_pH[14,5] 0.618 0.231 0.183 0.615 1.091
beta0_pH[15,5] 2.218 0.419 1.582 2.145 3.194
beta0_pH[16,5] -0.769 0.066 -0.897 -0.771 -0.635
beta1_pH[1,1] 2.602 0.500 1.816 2.527 3.834
beta1_pH[2,1] 1.575 0.305 1.021 1.563 2.187
beta1_pH[3,1] 1.800 0.415 1.154 1.754 2.866
beta1_pH[4,1] 1.943 0.450 1.196 1.887 3.099
beta1_pH[5,1] 2.376 0.507 1.562 2.317 3.602
beta1_pH[6,1] 2.303 0.857 1.148 2.099 4.378
beta1_pH[7,1] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,1] 2.684 1.033 1.228 2.477 5.349
beta1_pH[9,1] 2.126 0.631 1.243 2.009 3.748
beta1_pH[10,1] 2.217 0.400 1.527 2.170 3.136
beta1_pH[11,1] 3.303 0.556 2.367 3.218 4.574
beta1_pH[12,1] 2.876 0.286 2.331 2.878 3.466
beta1_pH[13,1] 3.427 0.233 2.981 3.420 3.893
beta1_pH[14,1] 2.791 0.337 2.193 2.779 3.478
beta1_pH[15,1] 2.373 0.279 1.879 2.357 2.970
beta1_pH[16,1] 4.137 0.317 3.532 4.124 4.767
beta1_pH[1,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[2,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[3,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[4,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[5,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[6,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[7,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[9,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[10,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[11,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[1,3] 3.452 1.637 1.316 2.966 7.254
beta1_pH[2,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[3,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[4,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[5,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[6,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[7,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,3] 4.780 2.231 2.219 4.054 10.980
beta1_pH[9,3] 8.444 2.664 4.210 8.026 14.147
beta1_pH[10,3] 6.441 2.548 2.854 5.924 12.187
beta1_pH[11,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[11,4] 2.114 1.012 0.430 1.992 5.000
beta1_pH[12,4] 4.027 1.098 2.089 3.986 6.506
beta1_pH[13,4] 2.539 0.813 1.164 2.470 4.307
beta1_pH[14,4] 2.758 0.734 1.550 2.679 4.451
beta1_pH[15,4] 1.169 0.991 0.113 0.829 3.659
beta1_pH[16,4] 3.888 0.916 2.381 3.769 6.024
beta1_pH[11,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[1,1] 1.222 0.845 0.323 0.915 3.375
beta2_pH[2,1] 1.631 0.953 0.335 1.404 3.422
beta2_pH[3,1] 1.369 0.855 0.194 1.191 2.942
beta2_pH[4,1] 1.088 0.784 0.205 0.809 2.910
beta2_pH[5,1] 2.118 1.915 0.253 1.569 7.649
beta2_pH[6,1] 0.965 1.627 0.103 0.352 6.458
beta2_pH[7,1] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,1] 0.827 1.645 0.111 0.288 6.818
beta2_pH[9,1] 1.234 1.849 0.152 0.494 7.310
beta2_pH[10,1] 1.324 1.592 0.274 0.753 6.786
beta2_pH[11,1] 0.429 0.257 0.178 0.370 1.077
beta2_pH[12,1] 1.597 0.814 0.636 1.345 3.695
beta2_pH[13,1] 1.031 0.432 0.542 0.923 2.210
beta2_pH[14,1] 1.111 0.655 0.436 0.906 3.121
beta2_pH[15,1] 1.080 0.650 0.377 0.900 3.115
beta2_pH[16,1] 1.104 0.412 0.630 1.009 2.236
beta2_pH[1,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[2,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[3,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[4,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[5,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[6,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[7,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[9,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[10,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[11,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[12,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[13,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[14,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[15,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[16,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[1,3] 2.178 3.808 0.111 0.491 14.803
beta2_pH[2,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[3,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[4,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[5,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[6,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[7,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,3] 2.422 1.938 0.584 1.651 7.914
beta2_pH[9,3] 0.188 0.165 0.081 0.159 0.429
beta2_pH[10,3] 0.381 0.486 0.109 0.265 1.410
beta2_pH[11,3] -8.215 5.925 -22.079 -6.962 -0.436
beta2_pH[12,3] -8.230 5.941 -22.106 -6.982 -0.411
beta2_pH[13,3] -8.322 5.993 -22.543 -7.061 -0.444
beta2_pH[14,3] -8.252 5.944 -22.410 -7.125 -0.326
beta2_pH[15,3] -8.311 5.901 -22.169 -7.073 -0.429
beta2_pH[16,3] -8.275 5.980 -22.091 -7.039 -0.360
beta2_pH[11,4] -4.301 3.741 -14.279 -3.099 -0.626
beta2_pH[12,4] -4.826 3.699 -15.175 -3.723 -0.893
beta2_pH[13,4] -4.056 3.672 -14.270 -2.750 -0.662
beta2_pH[14,4] -4.274 3.615 -14.255 -3.007 -0.803
beta2_pH[15,4] -4.451 3.665 -14.602 -3.207 -0.753
beta2_pH[16,4] -4.624 3.436 -14.546 -3.480 -1.236
beta2_pH[11,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[12,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[13,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[14,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[15,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[16,5] 0.000 0.000 0.000 0.000 0.000
beta3_pH[1,1] 34.980 1.209 32.760 34.950 37.817
beta3_pH[2,1] 32.646 1.438 29.967 32.582 35.714
beta3_pH[3,1] 33.049 1.597 30.241 32.963 36.862
beta3_pH[4,1] 32.668 1.857 29.649 32.457 37.064
beta3_pH[5,1] 27.572 1.181 25.854 27.437 30.656
beta3_pH[6,1] 32.898 3.640 25.322 32.959 40.827
beta3_pH[7,1] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,1] 34.927 3.838 28.626 34.702 42.906
beta3_pH[9,1] 32.475 2.981 27.643 31.935 39.549
beta3_pH[10,1] 32.767 1.377 30.425 32.686 35.838
beta3_pH[11,1] 31.216 1.121 29.234 31.166 33.577
beta3_pH[12,1] 30.793 0.477 29.776 30.813 31.678
beta3_pH[13,1] 31.340 0.539 30.339 31.334 32.418
beta3_pH[14,1] 32.600 0.776 31.255 32.553 34.269
beta3_pH[15,1] 30.898 0.719 29.478 30.912 32.300
beta3_pH[16,1] 30.084 0.330 29.420 30.089 30.725
beta3_pH[1,2] 34.451 5.466 25.443 34.400 43.479
beta3_pH[2,2] 34.592 5.432 25.527 34.639 43.505
beta3_pH[3,2] 30.399 7.234 18.640 30.420 42.361
beta3_pH[4,2] 30.502 7.114 18.680 30.343 42.262
beta3_pH[5,2] 30.333 7.207 18.552 30.274 42.342
beta3_pH[6,2] 30.325 7.166 18.625 30.568 42.199
beta3_pH[7,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[9,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[10,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[11,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[12,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[13,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[14,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[15,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[16,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[1,3] 39.793 2.520 33.542 40.163 43.693
beta3_pH[2,3] 34.647 5.494 25.611 34.594 43.500
beta3_pH[3,3] 30.476 7.323 18.644 30.352 42.443
beta3_pH[4,3] 30.432 7.236 18.679 30.248 42.331
beta3_pH[5,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[6,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[7,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,3] 42.097 0.807 40.977 41.853 43.828
beta3_pH[9,3] 35.956 1.745 34.527 35.285 41.422
beta3_pH[10,3] 37.914 2.457 34.710 37.312 43.324
beta3_pH[11,3] 38.517 2.614 34.188 38.449 42.799
beta3_pH[12,3] 38.450 2.600 34.223 38.441 42.750
beta3_pH[13,3] 38.512 2.619 34.221 38.464 42.778
beta3_pH[14,3] 38.529 2.578 34.242 38.495 42.734
beta3_pH[15,3] 38.572 2.587 34.261 38.608 42.827
beta3_pH[16,3] 38.571 2.574 34.262 38.563 42.768
beta3_pH[11,4] 40.751 1.319 38.274 40.660 43.375
beta3_pH[12,4] 42.237 0.669 40.430 42.380 43.101
beta3_pH[13,4] 40.962 0.944 39.241 40.914 42.985
beta3_pH[14,4] 40.245 0.735 38.854 40.224 41.677
beta3_pH[15,4] 41.553 1.672 37.531 41.804 43.823
beta3_pH[16,4] 41.733 0.460 40.830 41.777 42.685
beta3_pH[11,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[12,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[13,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[14,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[15,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[16,5] 25.000 0.000 25.000 25.000 25.000
beta4_pH[1,1] 0.262 0.021 0.221 0.262 0.303
beta4_pH[2,1] 0.211 0.005 0.201 0.211 0.221
beta4_pH[3,1] 0.173 0.007 0.160 0.173 0.187
beta4_pH[4,1] 0.294 0.011 0.272 0.293 0.316
beta4_pH[5,1] 0.422 0.272 0.092 0.383 0.947
beta4_pH[6,1] 0.775 0.380 0.332 0.734 1.436
beta4_pH[7,1] 0.848 0.331 0.360 0.797 1.644
beta4_pH[8,1] 0.153 0.069 0.039 0.147 0.308
beta4_pH[9,1] 0.215 0.087 0.058 0.210 0.390
beta4_pH[10,1] 0.001 0.001 0.000 0.001 0.003
beta4_pH[11,1] 0.446 0.012 0.423 0.447 0.470
beta4_pH[12,1] 0.643 0.022 0.599 0.642 0.688
beta4_pH[13,1] 0.630 0.014 0.604 0.631 0.656
beta4_pH[14,1] 0.670 0.031 0.610 0.669 0.730
beta4_pH[15,1] 0.797 0.021 0.758 0.797 0.839
beta4_pH[16,1] 1.210 0.042 1.131 1.210 1.295
beta4_pH[1,2] 0.297 0.095 0.126 0.290 0.494
beta4_pH[2,2] 0.515 0.023 0.471 0.515 0.558
beta4_pH[3,2] 0.446 0.012 0.423 0.446 0.470
beta4_pH[4,2] 0.656 0.038 0.584 0.654 0.733
beta4_pH[5,2] 0.731 0.880 0.102 0.527 2.707
beta4_pH[6,2] 0.907 0.870 0.188 0.740 2.639
beta4_pH[7,2] 0.897 0.420 0.366 0.815 1.932
beta4_pH[8,2] 0.849 0.203 0.488 0.843 1.251
beta4_pH[9,2] 1.120 1.063 0.297 0.861 3.511
beta4_pH[10,2] 0.457 0.050 0.362 0.456 0.558
beta4_pH[11,2] 1.268 0.050 1.174 1.267 1.366
beta4_pH[12,2] 0.608 0.114 0.400 0.603 0.841
beta4_pH[13,2] 0.955 0.044 0.870 0.954 1.045
beta4_pH[14,2] 0.936 0.063 0.815 0.937 1.063
beta4_pH[15,2] 1.011 0.026 0.961 1.011 1.064
beta4_pH[16,2] 1.369 0.071 1.236 1.367 1.509
beta4_pH[1,3] 0.240 0.076 0.099 0.235 0.399
beta4_pH[2,3] 0.049 0.017 0.018 0.049 0.085
beta4_pH[3,3] 0.052 0.010 0.033 0.052 0.072
beta4_pH[4,3] 0.167 0.031 0.107 0.167 0.227
beta4_pH[5,3] 0.760 0.862 0.115 0.556 2.702
beta4_pH[6,3] 0.895 0.787 0.204 0.725 2.506
beta4_pH[7,3] 0.893 0.418 0.346 0.818 1.909
beta4_pH[8,3] 0.385 0.185 0.093 0.361 0.818
beta4_pH[9,3] 0.928 0.825 0.163 0.722 2.764
beta4_pH[10,3] 0.056 0.038 0.004 0.050 0.142
beta4_pH[11,3] 0.621 0.263 0.258 0.574 1.255
beta4_pH[12,3] 0.738 0.267 0.361 0.693 1.375
beta4_pH[13,3] 0.793 0.374 0.296 0.725 1.684
beta4_pH[14,3] 0.872 0.224 0.514 0.848 1.382
beta4_pH[15,3] 0.794 0.202 0.463 0.770 1.263
beta4_pH[16,3] 0.828 0.370 0.333 0.764 1.767
beta4_pH[1,4] 0.637 0.376 0.198 0.561 1.518
beta4_pH[2,4] 0.803 0.833 0.112 0.589 2.734
beta4_pH[3,4] 0.860 0.927 0.114 0.605 3.089
beta4_pH[4,4] 0.733 0.464 0.209 0.632 1.936
beta4_pH[5,4] 0.734 0.749 0.104 0.557 2.500
beta4_pH[6,4] 0.917 0.791 0.202 0.736 2.708
beta4_pH[7,4] 0.895 0.416 0.341 0.820 1.896
beta4_pH[8,4] 0.765 0.827 0.145 0.591 2.413
beta4_pH[9,4] 0.774 0.928 0.117 0.553 2.790
beta4_pH[10,4] 1.426 2.872 0.007 0.383 11.036
beta4_pH[11,4] 0.427 0.019 0.390 0.426 0.465
beta4_pH[12,4] 0.535 0.044 0.448 0.533 0.624
beta4_pH[13,4] 0.836 0.033 0.770 0.836 0.900
beta4_pH[14,4] 0.876 0.060 0.766 0.874 0.995
beta4_pH[15,4] 0.617 0.018 0.581 0.617 0.654
beta4_pH[16,4] 0.664 0.043 0.583 0.663 0.751
beta4_pH[1,5] 0.640 0.368 0.201 0.558 1.544
beta4_pH[2,5] 0.874 1.033 0.118 0.593 3.410
beta4_pH[3,5] 0.868 1.031 0.114 0.597 3.110
beta4_pH[4,5] 0.743 0.617 0.201 0.621 2.025
beta4_pH[5,5] 0.742 0.733 0.108 0.563 2.553
beta4_pH[6,5] 0.921 0.811 0.206 0.743 2.913
beta4_pH[7,5] 0.888 0.433 0.341 0.815 1.926
beta4_pH[8,5] 0.728 0.593 0.132 0.589 2.134
beta4_pH[9,5] 0.809 1.000 0.127 0.585 2.909
beta4_pH[10,5] 1.508 3.136 0.006 0.390 11.828
beta4_pH[11,5] 0.308 0.082 0.154 0.305 0.473
beta4_pH[12,5] 0.874 0.287 0.466 0.819 1.595
beta4_pH[13,5] 2.557 0.352 1.906 2.550 3.276
beta4_pH[14,5] 0.856 0.202 0.521 0.839 1.307
beta4_pH[15,5] 0.927 0.076 0.801 0.920 1.098
beta4_pH[16,5] 0.951 0.382 0.388 0.886 1.879
beta0_pelagic[1] 2.267 0.036 2.201 2.266 2.341
beta0_pelagic[2] 1.538 0.007 1.524 1.538 1.551
beta0_pelagic[3] 0.138 0.273 -0.551 0.192 0.523
beta0_pelagic[4] 0.243 0.239 -0.317 0.253 0.676
beta0_pelagic[5] 1.438 0.154 1.148 1.436 1.746
beta0_pelagic[6] 1.356 0.059 1.246 1.355 1.473
beta0_pelagic[7] 1.468 0.039 1.395 1.467 1.544
beta0_pelagic[8] 1.671 0.043 1.589 1.670 1.758
beta0_pelagic[9] 2.773 0.063 2.650 2.772 2.900
beta0_pelagic[10] 2.604 0.031 2.547 2.603 2.666
beta0_pelagic[11] -0.122 0.184 -0.516 -0.109 0.228
beta0_pelagic[12] 1.800 0.035 1.734 1.799 1.873
beta0_pelagic[13] -1.002 0.529 -2.115 -0.973 -0.127
beta0_pelagic[14] -0.095 0.216 -0.608 -0.054 0.213
beta0_pelagic[15] -0.200 0.029 -0.280 -0.189 -0.174
beta0_pelagic[16] -0.821 0.238 -1.460 -0.789 -0.442
beta1_pelagic[1] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[2] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[3] 0.821 0.315 0.361 0.758 1.649
beta1_pelagic[4] 0.807 0.262 0.344 0.793 1.426
beta1_pelagic[5] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[6] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[7] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[8] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[9] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[10] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[11] 1.145 0.254 0.667 1.133 1.664
beta1_pelagic[12] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[13] 3.409 1.195 1.690 3.137 6.379
beta1_pelagic[14] 1.592 0.541 0.858 1.461 2.886
beta1_pelagic[15] 0.076 0.145 0.000 0.010 0.566
beta1_pelagic[16] 1.954 0.353 1.412 1.907 2.929
beta2_pelagic[1] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[2] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[3] 5.814 5.878 0.154 3.941 20.498
beta2_pelagic[4] 7.747 6.181 0.327 6.419 22.145
beta2_pelagic[5] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[6] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[7] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[8] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[9] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[10] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[11] 1.960 1.663 0.282 1.403 6.499
beta2_pelagic[12] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[13] 0.154 0.077 0.064 0.140 0.343
beta2_pelagic[14] 0.992 1.116 0.157 0.557 4.479
beta2_pelagic[15] 2.632 1.936 0.158 2.221 7.429
beta2_pelagic[16] 1.507 1.139 0.243 1.154 4.702
beta3_pelagic[1] 28.657 3.975 22.320 28.488 35.196
beta3_pelagic[2] 28.904 3.906 22.303 29.054 35.215
beta3_pelagic[3] 27.530 2.845 22.845 27.268 32.671
beta3_pelagic[4] 26.299 1.776 23.007 26.299 31.654
beta3_pelagic[5] 29.882 6.659 18.637 30.406 40.540
beta3_pelagic[6] 29.530 6.645 18.607 29.779 40.484
beta3_pelagic[7] 29.369 6.636 18.553 29.164 40.372
beta3_pelagic[8] 29.738 6.546 18.649 29.929 40.496
beta3_pelagic[9] 29.240 6.599 18.514 29.134 40.240
beta3_pelagic[10] 29.514 6.738 18.511 29.452 40.441
beta3_pelagic[11] 30.011 1.286 27.281 30.000 32.995
beta3_pelagic[12] 30.000 0.000 30.000 30.000 30.000
beta3_pelagic[13] 34.198 3.533 29.215 33.770 41.560
beta3_pelagic[14] 35.237 1.908 31.808 35.005 39.510
beta3_pelagic[15] 34.219 5.797 23.615 35.256 41.741
beta3_pelagic[16] 30.354 0.767 28.804 30.330 31.993
beta0_yellow[1] -0.487 0.027 -0.543 -0.487 -0.434
beta0_yellow[2] 0.283 0.143 -0.224 0.326 0.395
beta0_yellow[3] -0.297 0.080 -0.487 -0.285 -0.185
beta0_yellow[4] 0.623 0.255 -0.089 0.704 0.917
beta0_yellow[5] 0.119 0.332 -0.556 0.129 0.744
beta0_yellow[6] 0.836 0.062 0.715 0.835 0.959
beta0_yellow[7] 1.116 0.032 1.054 1.115 1.178
beta0_yellow[8] 1.050 0.029 0.992 1.050 1.107
beta0_yellow[9] 0.504 0.057 0.395 0.504 0.617
beta0_yellow[10] 0.181 0.028 0.128 0.181 0.237
beta0_yellow[11] -0.167 0.082 -0.378 -0.128 -0.093
beta0_yellow[12] -1.097 0.027 -1.149 -1.097 -1.045
beta0_yellow[13] -1.278 0.165 -1.676 -1.258 -1.048
beta0_yellow[14] -0.292 0.121 -0.550 -0.261 -0.167
beta0_yellow[15] -1.097 0.111 -1.418 -1.055 -1.018
beta0_yellow[16] -0.282 0.012 -0.305 -0.282 -0.259
beta1_yellow[1] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[2] 1.028 0.266 0.668 0.982 1.860
beta1_yellow[3] 0.570 0.153 0.277 0.571 0.877
beta1_yellow[4] 1.310 0.492 0.720 1.163 2.685
beta1_yellow[5] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[6] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[7] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[8] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[9] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[10] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[11] 0.258 0.349 0.001 0.063 1.146
beta1_yellow[12] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[13] 0.975 0.878 0.023 0.771 3.608
beta1_yellow[14] 0.618 0.901 0.013 0.360 3.499
beta1_yellow[15] 0.178 0.240 0.011 0.080 0.801
beta1_yellow[16] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[1] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[2] -5.383 5.693 -20.429 -3.151 -0.188
beta2_yellow[3] -7.893 5.947 -22.366 -6.454 -0.623
beta2_yellow[4] -4.323 5.510 -18.680 -1.551 -0.130
beta2_yellow[5] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[6] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[7] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[8] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[9] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[10] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[11] -6.489 5.314 -20.216 -4.967 -0.260
beta2_yellow[12] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[13] -6.112 5.436 -19.544 -4.508 -0.270
beta2_yellow[14] -7.338 5.619 -20.351 -6.174 -0.380
beta2_yellow[15] -7.359 5.466 -19.934 -6.410 -0.252
beta2_yellow[16] 0.000 0.000 0.000 0.000 0.000
beta3_yellow[1] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[2] 31.889 0.864 31.023 31.661 34.093
beta3_yellow[3] 32.803 1.480 30.301 32.667 36.104
beta3_yellow[4] 33.007 1.604 30.228 32.950 36.332
beta3_yellow[5] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[6] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[7] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[8] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[9] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[10] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[11] 35.159 3.706 30.382 34.177 43.273
beta3_yellow[12] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[13] 33.871 2.669 29.335 33.885 40.560
beta3_yellow[14] 34.284 3.672 29.207 33.793 43.104
beta3_yellow[15] 35.802 3.725 30.282 35.405 43.391
beta3_yellow[16] 35.000 0.000 35.000 35.000 35.000
beta0_black[1] -0.013 0.069 -0.153 -0.009 0.124
beta0_black[2] 1.852 0.031 1.794 1.852 1.916
beta0_black[3] 1.237 0.037 1.165 1.238 1.311
beta0_black[4] 2.412 0.046 2.322 2.413 2.502
beta0_black[5] 4.876 2.287 1.874 4.404 10.731
beta0_black[6] 4.713 2.067 2.310 4.197 10.131
beta0_black[7] 3.463 1.950 1.344 2.947 9.062
beta0_black[8] 0.703 0.171 0.382 0.698 1.057
beta0_black[9] 2.641 0.228 2.219 2.636 3.092
beta0_black[10] 1.276 0.049 1.184 1.276 1.377
beta0_black[11] 3.388 0.024 3.342 3.388 3.434
beta0_black[12] 4.925 0.099 4.740 4.924 5.127
beta0_black[13] 0.243 0.190 0.009 0.200 0.693
beta0_black[14] 3.005 0.041 2.927 3.005 3.086
beta0_black[15] 1.242 0.026 1.190 1.242 1.294
beta0_black[16] 4.261 0.057 4.151 4.261 4.376
beta2_black[1] 0.000 0.000 0.000 0.000 0.000
beta2_black[2] 0.000 0.000 0.000 0.000 0.000
beta2_black[3] 0.000 0.000 0.000 0.000 0.000
beta2_black[4] 0.000 0.000 0.000 0.000 0.000
beta2_black[5] 0.000 0.000 0.000 0.000 0.000
beta2_black[6] 0.000 0.000 0.000 0.000 0.000
beta2_black[7] 0.000 0.000 0.000 0.000 0.000
beta2_black[8] 0.000 0.000 0.000 0.000 0.000
beta2_black[9] 0.000 0.000 0.000 0.000 0.000
beta2_black[10] 0.000 0.000 0.000 0.000 0.000
beta2_black[11] 0.000 0.000 0.000 0.000 0.000
beta2_black[12] 0.000 0.000 0.000 0.000 0.000
beta2_black[13] -2.519 1.638 -6.539 -2.097 -0.425
beta2_black[14] 0.000 0.000 0.000 0.000 0.000
beta2_black[15] 0.000 0.000 0.000 0.000 0.000
beta2_black[16] 0.000 0.000 0.000 0.000 0.000
beta3_black[1] 25.000 0.000 25.000 25.000 25.000
beta3_black[2] 25.000 0.000 25.000 25.000 25.000
beta3_black[3] 25.000 0.000 25.000 25.000 25.000
beta3_black[4] 25.000 0.000 25.000 25.000 25.000
beta3_black[5] 25.000 0.000 25.000 25.000 25.000
beta3_black[6] 25.000 0.000 25.000 25.000 25.000
beta3_black[7] 25.000 0.000 25.000 25.000 25.000
beta3_black[8] 25.000 0.000 25.000 25.000 25.000
beta3_black[9] 25.000 0.000 25.000 25.000 25.000
beta3_black[10] 25.000 0.000 25.000 25.000 25.000
beta3_black[11] 25.000 0.000 25.000 25.000 25.000
beta3_black[12] 25.000 0.000 25.000 25.000 25.000
beta3_black[13] 38.613 1.278 35.877 38.859 40.253
beta3_black[14] 25.000 0.000 25.000 25.000 25.000
beta3_black[15] 25.000 0.000 25.000 25.000 25.000
beta3_black[16] 25.000 0.000 25.000 25.000 25.000
beta4_black[1] -0.296 0.085 -0.462 -0.298 -0.125
beta4_black[2] 0.257 0.052 0.155 0.259 0.358
beta4_black[3] -0.736 0.059 -0.851 -0.736 -0.620
beta4_black[4] 0.543 0.119 0.312 0.545 0.774
beta4_black[5] 0.730 1.748 -1.011 0.451 3.773
beta4_black[6] 0.772 1.486 -0.728 0.489 3.799
beta4_black[7] 0.716 1.597 -0.863 0.448 3.677
beta4_black[8] 0.079 0.250 -0.450 0.083 0.556
beta4_black[9] 0.925 0.764 -0.131 0.772 2.883
beta4_black[10] 0.084 0.066 -0.047 0.084 0.211
beta4_black[11] -0.704 0.050 -0.799 -0.703 -0.606
beta4_black[12] 0.896 0.234 0.462 0.889 1.387
beta4_black[13] -1.174 0.048 -1.265 -1.174 -1.082
beta4_black[14] -0.063 0.113 -0.276 -0.064 0.157
beta4_black[15] -0.481 0.036 -0.552 -0.481 -0.414
beta4_black[16] -0.597 0.108 -0.805 -0.597 -0.374
beta0_dsr[11] -0.338 0.009 -0.356 -0.338 -0.320
beta0_dsr[12] 1.421 0.043 1.338 1.420 1.505
beta0_dsr[13] 0.142 0.021 0.100 0.142 0.185
beta0_dsr[14] -0.128 0.017 -0.162 -0.129 -0.095
beta0_dsr[15] 0.197 0.014 0.168 0.197 0.224
beta0_dsr[16] 0.208 0.015 0.179 0.208 0.237
beta1_dsr[11] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[12] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[13] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[14] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[15] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[16] 0.000 0.000 0.000 0.000 0.000
beta4_dsr[11] 0.549 0.035 0.483 0.548 0.616
beta4_dsr[12] 0.789 0.232 0.350 0.787 1.241
beta4_dsr[13] -0.079 0.040 -0.158 -0.079 0.001
beta4_dsr[14] -0.214 0.073 -0.358 -0.214 -0.071
beta4_dsr[15] 0.225 0.022 0.183 0.226 0.267
beta4_dsr[16] -0.100 0.044 -0.187 -0.101 -0.014
beta0_slope[11] -2.198 0.021 -2.239 -2.198 -2.158
beta0_slope[12] -5.157 0.175 -5.492 -5.159 -4.802
beta0_slope[13] -1.473 0.034 -1.538 -1.473 -1.407
beta0_slope[14] -2.685 0.048 -2.782 -2.683 -2.594
beta0_slope[15] -1.589 0.022 -1.633 -1.590 -1.545
beta0_slope[16] -3.014 0.046 -3.106 -3.014 -2.923
beta1_slope[11] 0.000 0.000 0.000 0.000 0.000
beta1_slope[12] 0.000 0.000 0.000 0.000 0.000
beta1_slope[13] 0.000 0.000 0.000 0.000 0.000
beta1_slope[14] 0.000 0.000 0.000 0.000 0.000
beta1_slope[15] 0.000 0.000 0.000 0.000 0.000
beta1_slope[16] 0.000 0.000 0.000 0.000 0.000
beta4_slope[11] -0.345 0.076 -0.496 -0.346 -0.198
beta4_slope[12] -4.235 0.798 -5.928 -4.204 -2.750
beta4_slope[13] 0.305 0.053 0.201 0.304 0.411
beta4_slope[14] -0.186 0.154 -0.489 -0.184 0.117
beta4_slope[15] -0.364 0.034 -0.431 -0.364 -0.297
beta4_slope[16] 0.218 0.104 0.012 0.217 0.415
beta5_slope[11] 3.078 0.108 2.876 3.074 3.295
beta5_slope[12] 1.571 0.410 0.810 1.554 2.416
beta5_slope[13] 0.865 0.079 0.717 0.864 1.021
beta5_slope[14] 3.776 0.357 3.151 3.753 4.554
beta5_slope[15] 1.384 0.073 1.247 1.383 1.530
beta5_slope[16] 2.724 0.231 2.288 2.718 3.199
tau_comp[1] 3.850 0.305 3.293 3.838 4.491
tau_comp[2] 2.881 0.256 2.391 2.875 3.399
tau_comp[3] 1.653 0.133 1.405 1.648 1.925
tau_comp[4] 1.344 0.211 0.980 1.328 1.799
tau_comp[5] 1.210 0.180 0.897 1.199 1.592
sigma_H[1] 0.206 0.049 0.121 0.203 0.315
sigma_H[2] 0.175 0.028 0.126 0.173 0.235
sigma_H[3] 0.188 0.040 0.117 0.184 0.275
sigma_H[4] 0.464 0.081 0.331 0.455 0.653
sigma_H[5] 0.963 0.201 0.595 0.955 1.397
sigma_H[6] 0.523 0.191 0.143 0.514 0.917
sigma_H[7] 0.310 0.062 0.211 0.301 0.454
sigma_H[8] 0.421 0.088 0.283 0.411 0.611
sigma_H[9] 0.513 0.116 0.330 0.497 0.784
sigma_H[10] 0.209 0.038 0.144 0.206 0.294
sigma_H[11] 0.266 0.043 0.194 0.262 0.361
sigma_H[12] 0.382 0.148 0.199 0.336 0.725
sigma_H[13] 0.221 0.036 0.160 0.218 0.300
sigma_H[14] 0.495 0.085 0.342 0.490 0.676
sigma_H[15] 0.228 0.038 0.164 0.225 0.310
sigma_H[16] 0.214 0.041 0.147 0.211 0.304
lambda_H[1] 2.201 3.270 0.073 1.074 11.239
lambda_H[2] 5.286 5.007 0.517 3.715 18.836
lambda_H[3] 5.577 8.219 0.266 2.685 28.511
lambda_H[4] 0.006 0.004 0.001 0.005 0.018
lambda_H[5] 3.065 8.311 0.026 0.703 24.517
lambda_H[6] 5.676 11.971 0.009 0.242 39.860
lambda_H[7] 0.012 0.009 0.002 0.010 0.033
lambda_H[8] 8.020 10.507 0.099 4.336 37.812
lambda_H[9] 0.015 0.010 0.003 0.013 0.040
lambda_H[10] 0.225 0.320 0.025 0.159 0.825
lambda_H[11] 0.236 0.331 0.010 0.113 1.077
lambda_H[12] 3.765 5.368 0.164 2.139 16.962
lambda_H[13] 3.150 3.046 0.266 2.254 11.003
lambda_H[14] 3.392 4.179 0.224 2.124 13.691
lambda_H[15] 0.018 0.019 0.003 0.014 0.062
lambda_H[16] 0.563 0.973 0.037 0.299 2.616
mu_lambda_H[1] 4.045 1.855 1.027 3.869 8.081
mu_lambda_H[2] 3.639 1.955 0.530 3.454 7.797
mu_lambda_H[3] 3.369 1.811 0.715 3.092 7.483
sigma_lambda_H[1] 8.077 4.277 1.628 7.351 17.808
sigma_lambda_H[2] 7.976 4.718 0.877 7.328 18.250
sigma_lambda_H[3] 6.166 3.977 0.961 5.287 16.377
beta_H[1,1] 6.471 1.316 3.268 6.680 8.391
beta_H[2,1] 9.764 0.624 8.450 9.799 10.944
beta_H[3,1] 7.817 0.813 5.901 7.916 9.158
beta_H[4,1] 10.671 7.597 -4.859 10.666 25.897
beta_H[5,1] -0.250 2.545 -5.539 -0.120 4.036
beta_H[6,1] 2.958 4.460 -7.716 4.535 8.068
beta_H[7,1] 1.086 5.858 -11.764 1.529 11.651
beta_H[8,1] 1.260 3.417 -2.454 1.175 3.405
beta_H[9,1] 13.646 5.595 2.764 13.585 24.958
beta_H[10,1] 6.799 1.807 2.951 6.914 10.234
beta_H[11,1] 4.911 3.546 -3.139 5.543 9.893
beta_H[12,1] 2.609 1.118 0.572 2.525 5.020
beta_H[13,1] 9.089 0.936 7.215 9.181 10.611
beta_H[14,1] 2.133 1.021 0.045 2.143 4.203
beta_H[15,1] -7.124 3.393 -13.543 -7.181 -0.247
beta_H[16,1] 3.749 2.966 -0.826 3.276 10.650
beta_H[1,2] 7.791 0.272 7.249 7.791 8.311
beta_H[2,2] 9.878 0.162 9.555 9.880 10.199
beta_H[3,2] 8.798 0.202 8.417 8.795 9.210
beta_H[4,2] 3.042 1.439 0.277 3.045 5.953
beta_H[5,2] 1.747 1.002 -0.210 1.758 3.664
beta_H[6,2] 5.988 1.091 3.562 6.132 7.738
beta_H[7,2] 2.432 1.089 0.547 2.331 4.843
beta_H[8,2] 3.000 1.016 1.358 3.110 4.189
beta_H[9,2] 3.191 1.079 1.096 3.172 5.338
beta_H[10,2] 8.043 0.360 7.299 8.057 8.732
beta_H[11,2] 9.818 0.626 8.880 9.717 11.219
beta_H[12,2] 3.937 0.372 3.231 3.923 4.700
beta_H[13,2] 9.144 0.244 8.700 9.137 9.621
beta_H[14,2] 4.010 0.337 3.333 4.011 4.683
beta_H[15,2] 11.560 0.622 10.281 11.576 12.743
beta_H[16,2] 4.499 0.817 3.002 4.482 6.169
beta_H[1,3] 8.448 0.260 7.981 8.434 8.972
beta_H[2,3] 9.854 0.147 9.568 9.853 10.162
beta_H[3,3] 9.465 0.169 9.153 9.459 9.813
beta_H[4,3] -2.189 0.856 -3.869 -2.171 -0.569
beta_H[5,3] 3.878 0.614 2.641 3.872 5.051
beta_H[6,3] 8.364 1.203 6.624 8.056 10.849
beta_H[7,3] -2.576 0.637 -3.868 -2.563 -1.380
beta_H[8,3] 5.264 0.487 4.668 5.204 6.321
beta_H[9,3] -2.732 0.705 -4.098 -2.727 -1.370
beta_H[10,3] 8.547 0.273 8.028 8.547 9.085
beta_H[11,3] 8.501 0.289 7.869 8.528 9.012
beta_H[12,3] 5.263 0.325 4.488 5.296 5.799
beta_H[13,3] 8.855 0.180 8.483 8.862 9.190
beta_H[14,3] 5.770 0.267 5.159 5.789 6.250
beta_H[15,3] 10.268 0.292 9.718 10.263 10.850
beta_H[16,3] 6.107 0.565 4.970 6.125 7.101
beta_H[1,4] 8.103 0.184 7.717 8.119 8.428
beta_H[2,4] 9.850 0.130 9.585 9.856 10.087
beta_H[3,4] 9.946 0.155 9.613 9.956 10.217
beta_H[4,4] 11.808 0.473 10.879 11.808 12.722
beta_H[5,4] 5.562 0.731 4.355 5.493 7.209
beta_H[6,4] 7.137 0.909 5.233 7.336 8.460
beta_H[7,4] 8.555 0.350 7.841 8.561 9.238
beta_H[8,4] 6.719 0.252 6.228 6.731 7.139
beta_H[9,4] 7.436 0.452 6.572 7.435 8.320
beta_H[10,4] 7.523 0.218 7.097 7.520 7.966
beta_H[11,4] 9.496 0.194 9.132 9.496 9.875
beta_H[12,4] 7.251 0.198 6.852 7.243 7.672
beta_H[13,4] 9.053 0.144 8.758 9.057 9.325
beta_H[14,4] 7.818 0.212 7.414 7.812 8.248
beta_H[15,4] 9.472 0.221 9.048 9.472 9.908
beta_H[16,4] 9.440 0.233 9.002 9.433 9.914
beta_H[1,5] 8.906 0.154 8.584 8.914 9.200
beta_H[2,5] 10.679 0.101 10.485 10.676 10.885
beta_H[3,5] 10.760 0.163 10.469 10.753 11.091
beta_H[4,5] 8.358 0.498 7.397 8.344 9.389
beta_H[5,5] 5.186 0.604 3.728 5.257 6.160
beta_H[6,5] 8.968 0.647 7.988 8.856 10.351
beta_H[7,5] 6.567 0.347 5.911 6.560 7.282
beta_H[8,5] 8.249 0.206 7.898 8.234 8.677
beta_H[9,5] 8.194 0.460 7.286 8.187 9.100
beta_H[10,5] 10.042 0.220 9.606 10.041 10.485
beta_H[11,5] 11.586 0.225 11.147 11.581 12.051
beta_H[12,5] 8.603 0.187 8.243 8.604 8.982
beta_H[13,5] 10.133 0.134 9.881 10.128 10.408
beta_H[14,5] 9.241 0.216 8.826 9.234 9.703
beta_H[15,5] 11.299 0.230 10.840 11.302 11.740
beta_H[16,5] 9.873 0.180 9.511 9.877 10.214
beta_H[1,6] 10.122 0.246 9.740 10.085 10.741
beta_H[2,6] 11.470 0.113 11.244 11.471 11.688
beta_H[3,6] 10.620 0.155 10.293 10.628 10.904
beta_H[4,6] 12.905 0.879 11.146 12.928 14.575
beta_H[5,6] 5.672 0.653 4.551 5.622 7.027
beta_H[6,6] 8.716 0.576 7.372 8.750 9.823
beta_H[7,6] 10.577 0.600 9.397 10.580 11.733
beta_H[8,6] 9.535 0.261 9.051 9.548 9.968
beta_H[9,6] 8.301 0.759 6.880 8.275 9.861
beta_H[10,6] 9.141 0.334 8.437 9.161 9.746
beta_H[11,6] 10.695 0.365 9.884 10.734 11.306
beta_H[12,6] 9.381 0.241 8.921 9.372 9.898
beta_H[13,6] 11.064 0.156 10.766 11.063 11.377
beta_H[14,6] 9.917 0.277 9.354 9.920 10.436
beta_H[15,6] 10.530 0.417 9.733 10.517 11.363
beta_H[16,6] 10.587 0.244 10.074 10.601 11.049
beta_H[1,7] 10.148 1.304 6.672 10.468 11.950
beta_H[2,7] 12.120 0.503 10.999 12.144 13.052
beta_H[3,7] 10.277 0.677 8.750 10.337 11.420
beta_H[4,7] 1.489 4.511 -7.107 1.498 10.519
beta_H[5,7] 5.821 2.153 1.773 5.926 9.776
beta_H[6,7] 7.156 3.191 -1.026 8.218 11.091
beta_H[7,7] 7.525 2.956 1.800 7.512 13.367
beta_H[8,7] 10.893 0.936 9.159 10.885 12.715
beta_H[9,7] 4.916 3.901 -2.913 4.953 12.644
beta_H[10,7] 10.028 1.559 7.360 9.891 13.445
beta_H[11,7] 11.335 1.812 8.257 11.120 15.366
beta_H[12,7] 9.777 0.979 7.523 9.864 11.399
beta_H[13,7] 11.680 0.708 10.092 11.735 12.886
beta_H[14,7] 10.642 0.882 8.840 10.680 12.271
beta_H[15,7] 14.058 2.185 9.683 14.066 18.292
beta_H[16,7] 12.829 1.323 10.610 12.682 15.774
beta0_H[1] 8.951 18.434 -25.516 8.407 43.195
beta0_H[2] 10.116 7.793 -6.164 10.189 26.029
beta0_H[3] 9.998 9.737 -8.899 9.805 30.476
beta0_H[4] 7.290 181.459 -349.131 3.599 379.088
beta0_H[5] 4.706 30.307 -53.919 4.147 60.045
beta0_H[6] 5.804 52.110 -111.226 7.663 113.106
beta0_H[7] 4.707 140.994 -269.620 4.959 287.264
beta0_H[8] 6.204 20.941 -19.273 6.315 30.632
beta0_H[9] 6.059 123.322 -246.866 6.918 243.720
beta0_H[10] 7.257 36.913 -67.003 7.423 84.302
beta0_H[11] 9.646 52.661 -105.220 10.214 119.353
beta0_H[12] 6.976 12.450 -16.430 6.823 31.393
beta0_H[13] 9.564 10.954 -12.829 9.626 30.754
beta0_H[14] 6.787 11.191 -16.623 7.013 28.870
beta0_H[15] 7.521 118.431 -231.470 6.867 245.420
beta0_H[16] 8.726 29.199 -49.025 7.783 72.520